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Abstract 


Several significant accomplishments were made during the present reporting period. 

• Initial simulations to understand the applicability of the MODIS 
1380 run band for removing the effects of stratospheric aerosols and 
thin cirrus clouds were completed using a model for an aged vol- 
canic aerosol. The results suggest that very simple procedures re- 
quiring no a priori knowledge of the optical properties of the strato- 
spheric aerosol may be as effective as complex procedures requiring 
full knowledge of the aerosol properties, except the concentration 
which is estimated from the reflectance at 1380 run. The limitations 
of this conclusion will be examined in the next reporting period. ■ 

• The lookup tables employed in the implementation of the atmo- 
spheric correction algorithm have been modified in several ways in- 
tended to improve the accuracy and/or speed of processing. These 
have been delivered to R. Evans for implementation into the MODIS 
prototype processing algorithm for testing. 

• A method was developed for removal of the effects of the O 2 “A" 
absorption band from SeaWiFS band 7 (745-7S5 nm). This is im- 
portant in that SeaWiFS imagery will be used as a test data set for 
the MODIS atmospheric correction algorithm over the oceans, 

• Construction of a radiometer, and associated deployment boom, for 
studying the spectral reflectance of oceanic whitecaps at sea was 
completed. The system w T as successfully tested on a cruise off Hawaii 
on which whitecaps were plentiful during October-November. This 
data set is now under analysis. 
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1. Atmospheric Correction Algorithm Development 

a. Task Objectives: 

During CY 1994 there are five objectives under this task: 

(i) Investigate the effects of stratospheric aerosol on the proposed correction algorithm, and 
investigate the possible use of the 1380 run MODIS band to remove the stratospheric aerosol 
perturbation. 

(ii) Investigate the effect of vertical structure in aerosol concentration and type on the behavior 
of the proposed correction algorithm. 

(iii) Investigate the effects of polarization on the accuracy of the algorithm. 

(iv) Improve the accuracy and speed of the existing algorithm. 

(v) Investigate removal of the O 2 “A 1 ’ absorption band at 762 run from the 745-785 nm 
SeaWiFS band so the latter can be used in atmospheric correction of SeaWiFS. The importance 
of this to MODIS is that SeaWiFS data will be used extensively to test and improve the MODIS 
algorithm. Thus, it is essential for the MODIS algorithm development and testing, that the O 2 
absorption be adequately dealt with for SeaWiFS. 

b. Work Accomplished: 

(i) We have completed the addition of a third layer to our successive order of scattering 
radiative transfer code to include stratospheric aerosols. A literature search was carried out and 
three stratospheric aerosol models were selected representing a background stratospheric aerosol, 
a fresh volcanic aerosol, and an aged volcanic aerosol. Computations of the optical properties of 
the aerosol were carried out and initial radiative transfer simulations have been performed. The 
third layer in the model contains the stratospheric aerosol. To simulate the radiance observed in 
the 1380 nm band we use a single layer in cult code and omit the Fresnel-reflecting sea surface. 
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This in effect assumes that all of the photons that penetrate through the stratospheric aerosol are 
immediately absorbed by water vapor in the free troposphere. 

In the absence of stratospheric aerosol, the total reflectance p>(A), at a wavelength A, measured 
by MODIS can be decomposed as follows: 

p t (A) = p r (A) + + Pra(^) + (l) 

where p r (A) is the radiance resulting from multiple scattering by air molecules (Rayleigh scattering) 
in the absence of aerosols, p a (A) is the radiance resulting from multiple scattering by aerosols in the 
absence of the air, and p ra (A) is the interaction term between molecular and aerosol scattering. 1 In 
this equation, t is the diffuse transmittance of the atmosphere along the viewing direction specified 
by 9 Vi the angle between the normal to the sea surface and the sensor.* Radiance arising from 
specular reflection of direct sunlight from the sea surface (sun glitter) has been ignored because 
SeaWiFS can be tilted away from the glitter pattern. The influence of whitecaps has also been 
ignored under the assumption that their contribution can be removed from an estimate of the 
surface wind speed. 3 The goal of the atmospheric correction is the retrieval of p w from p t • 

The stratospheric aerosol contributes to the reflectance at all oi the MODIS ocean bands. 
Thus, in the presence of the stratospheric aerosol layer the reflectance will be increased by an 
amount 6p\*\ i.e. , 

P ( / ) (A) = p t (A) + ^i J) (A), 

where p[ ^ is the reflectance of the entire ocean- atmosphere system in the presence of stratospheric 
aerosol. As much of this contribution as possible should be removed from the visible and NIR bands 
before applying the atmospheric correction. Thus, the goal is to be able to remove fp^(A) from 
p^(A) before applying the atmospheric correction algorithm, which was developed for application 
to p t ( A). 

We have carried out simulations of radiative transfer with our three-layer model to examine 
several possibilities for utilizing the 1380 nm band for correction for stratospheric aerosols. Al- 
though, the analysis of these computations is incomplete, we will provide a brief report of the 
results thus far. 
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It is assumed that the reflectance at 1380 n in is totally due to the stratospheric aerosol. Also, 
as usual it is assumed that the aerosol properties in the lower layer are completely unknown. The 
atmospheric correction algorithm, as presently implemented, is then operated in several ways using 
the three-layer computations as pseudodata. 

1. The “measured” reflectances at 443, 765 and 865 are used in the al- 
gorithm as usual, i.e., no attention is paid to the fact that a strato- 
spheric aerosol may be present [p^( A) assumed to be p : (A)], and 
the error in the atmospheric correction at 443 nm is determined. 

2. The stratospheric aerosol is incorporated into the algorithm by sim- 
ply subtracting the reflectance at 1380 nm from those at 44 3, 765, 

and 865, i.e., p t (A) = /> ( ( j) (A) - J - 1 ( 1 3S0 ) . These are then inserted 

into the correction algorithm and determining the error in the cor- 
rection at 443 nm is determined. 


3. It is assumed that the spectral variation of the optical thickness of 
the stratospheric aerosol is known (either from an instrument, like 
SAGE or from measurements from the surface). The reflectance 
at 1380 nm (due entirely to the stratospheric aerosol) is scaled by 
the ratio of the stratospheric optical depth at the given A. r,(A), to 
that at (or in the case of surface measurements, near) 13S0 nm, and 
subtracted from the measured reflectances at the other wavelengths, 


i.e., 


p t (X) = p\ 9 \x) 


' s ^ ' * ~ 3 • ( 1 
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The Pt(A) are then inserted into the correction algorithm and the 
error in the correction at 443 nm determined. 


4. It is assumed that accurate measurements or predictions of the other 
optical properties of the stratospheric aerosol, the spectral scatter- 
ing phase function and single scattering albedo, along with the spec- 
tral variation of the optical depth are available for the stratospheric 
aerosol. The reflectance at 1380 nm is then scaled, by the ratio of 
the single-scattered stratospheric aerosol reflectances at A to that 
at 1380 nm, and subtracted from the reflectances in the visible and 
NIR, i.e., 


Pt{ A) = P ( , 0 (A) - 


t^(A)T,(A)p < (fl v ,<g v ;g 0 ,<fto; A) r (j) M3S0) 


where 


p» (*V 7 4>v 7 *0 7 


<£o;A) = J>,(0-,A) + (r(0 B ) + r(0 o ))P,(0+,A), 


cos 6± = ± cos 8 0 cos 9 V - sin 6q sin 9 V cos(<p v - <po)i 
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and r(a) is the Fresnel reflectance of the interface for an incident 
angle a. The parameters r J (A),^ J (A), and P,(a, A) are, respectively, 
the stratospheric aerosol optical thickness, single scattering albedo, 
and scattering phase function for a scattering angle a. The angles 0 o 
and 4>q are, respectively, the zenith and azimuth angles of a vector 
from the point on the sea surface under examination (pixel) to the 
sun, and likewise, 8 V and 4> v are the zenith and azimuth angles of 
a vector from the pixel to the sensor. These are measured with 
respect to the upward normal so 8 V and $q are both less than 90°. 
The resulting values of p* ( A ) are then inserted into the correction 
algortihm and the error at 443 nm is determined. Thus procedure is 
based on the assumption that the stratospheric aerosol enhancement 
of p t is all due to single scattering. 

5. As in (4) it is assumed that all of the optical properties of the aerosol 
are known except the concentration. A one-layer multiple scattering 
code (with a totally absorbing lower surface to represent the tropo- 
sphere) is used to determine r 3 (13S0) from p. 3 ' ( 1380). This deter- 
mines all of the optical properties cf the stratospheric aerosol. These 
properties are inserted into a one-layer multiple scattering code (with 
a Fresnel-reflecting sea surface as the lower boundary) to compute 

£/>^(A), which is subtracted from the measured reflectances p l : 5 ^(A) 
to provide p t ( A). The resulting values of p : .[ A) are then inserted 
into the correction algortihm and the error at 443 nm is determined. 
This procedure is based on the assumption that there is no radiative 
interaction between the stratospheric aerosol layer and the otner two 
layers. 

6. Except for the step in which £///'( A) is removed from p ; (A), this is 
identical to procedure 5. Once all of the optical properties of the 
stratospheric aerosol are known, they are inserted into a two-layer 
multiple scattering code (as opposed to a one-layer code in procedure 
5 above) with a Fresnel-reflecting sea surface as the lower boundary. 
The top layer consists of the stratospheric aerosol and the lower 
layer has only Rayleigh scattering. This incorporates the Rayleigh- 
stratospheric aerosol interaction explicitly (albeit approximately be- 
cause of the absence of the tropospheric aerosol), leaving only the 
Rayleigh-tropospheric aerosol and stratospheric-tropospheric aerosol 
interactions unaddressed. After subtracting the result of this compu- 
tation from p\*\ A), the result is inserted into the standard correction 
algorithm in which allowanceis made for the fact that p r has already 
been removed with the stratospheric component. This approach is 
possible because the properties of the Rayleigh scattering layer are 
completely known. 


We now present the results obtained using these procedures for cases with an aged volcanic 
aerosol in the stratosphere with r , = 0.05 and 0.15, and a Maritime aerosol (relative humidity 
80%) with r a = 0.15 in the marine boundary layer. In the absence of stratospheric aerosol, the 
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algorithm yields an almost perfect atmospheric correction. Figures 1 and 2 provide the error in the 
recovered value of tp w (443) for each of the procedures above. On the figures the procedure number 
is given with the description of the symbols. The goal is that the error approach that which would 
be obtained were the stratospheric aerosol absent, which is provided by the very heavy black line. 
In the model description on the figures, the symbols U 3L,” ;i 2L,” “R,” “S,” and U (R+ S)” stand for 
radiances computed from a three-layer model, a two-layer model, a one-layer model with Rayleigh 
scattering only, a one-layer model with stratospheric aerosol scattering only, and a two-layer model 
with stratospheric aerosols on the top and molecular scattering in the bottom. 

To no surprise, the most complex method of dealing with the stratospheric aerosol (No. 6 
above) yields the best overall correction; however, the second best methods (Nos. 1 and 2) are 
the simplest and require absolutely no a propn knowledge of the stratospheric aerosol’s optical 
properties. This is similar the Gordon and Castano 4 conclusion regarding the influence of the El 
Chichon aerosol on CZCS atmospheric correction, i.e.. that the presence of the aerosol had little 
effect on atmospheric correction. 

(ii) A Monte Carlo code employing a fifty-layer atmosphere has been developed to carry out 
this investigation. To realistically treat the aerosol in the code, we divided the atmosphere into 
four broad regions: (1) the marine boundary layer from the surface to 2 km, where the aerosol 
concentration is independent of altitude; (2) the free troposphere, where the aerosol concentrations 
varies in proportion to exp [—z/h], where z is the altitude (2-12 km) and h (the scale height) is 2 km; 
(3) the background stratosphere (12-30 km), w T here the aerosol concentration is also exponential 
with a scale height of 5 km; and (4) a volcanic region (20-25 km) within the stratosphere w r hich 
can contain a uniformly mixed volcanic aerosol. The optical properties of each of the four regions 
can be characterized by individual aerosol models, and any of the regions can be free of aerosols if 
desired. 


QWNL * 
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Aged Volcanic-- Viewing at Center, r, = 0.15, r #t = 0.05 
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Figure 1. Error in retrieved water-leaving reflectance as a function of the solar zenith angle for 
several methods of dealing with the stratospheric aerosol (1-6). r a = 0.15, r j£ = 0.05. Top: scan 
center, bottom: scan edge. 
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Figure 2. Error in retrieved water-leaving reflectance as a function of the solar zenith angle for 
several methods of dealing with the stratospheric aerosol (1-6). r a = 0.15, r jt = 0.15. Top: scan 
center, bottom: scan edge. 
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The models for the two lower regions are taken from Shettle and Fenn . 5 Briefly, based on 
size distribution and composition measurements they developed two models with log-normal size 
distributions called the Tropospheric (to model the aerosol in the free troposphere few large 
particles) and the Oceanic (to model the aerosol produced by sea spray few small particles). 
They combined these to form- the Maritime model to represent the aerosol in the marine boundary 
layer. In terms of total aerosol number per unit volume t lie Maritime model consists of 99% 
Tropospheric and 1 % Oceanic. Gordon and Wang added a Coastal model (99.5% Tropospheric 
and 0.5% Oceanic) to provide a description of the aerosol that might be more representative of the 
boundary layer near the coast. 

For the two upper regions we use a model for the background stratosphere from the WMO 7 
and for the volcanic aerosol from King et al. s Both models assume a 75% solution of H 2 SO 4 . We 
have also included a volcanic ash model * to represent fresh volcanic aerosol. The size distributions 
in these regions are modified gamma distributions. 

The code contains provision for a wind-roughened sea surface with a surface slope distribution 
governed by the Cox an Munk 9 probability density function. 

Preliminary computations (using the code) relating to the influence of vertical structure of 
the aerosol on atmospheric correction (the algorithm assumes that the aerosol is all in the surface 
boundary layer) have been carried out. These show an increased error in the retrievals in the 
presence of vertical structure; however, the error is not excessive. We have not progressed further 
with this aspect of the study as we decided to complete the addition of polarization to the code 
first (the next item). 

(iii) We have added polarization to our fifty-layer Monte Carlo code. Initial exercising of the 
code has been for an aerosol-free atmsophere. In this case the results agree with our successive order 
of scattering code^to within ^ 0.1% for an atmosphere bounded by a fiat Fresnel-reflecting ocean. 
For a sea surface roughened by the wind, differences in the two codes of the order of — 0.5% persist.. 
Considerable effort on this part of the project during the last three months has been devoted to 
tracking down the source of this difference. 
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(iv) For the implementation of our correction algorithm, 6 extensive lookup tables relating 
Pt ~~ Pr — tp w = p a + p ra to the single-scattered aerosol radiance, p as , are required for each aerosol 
model we employ. Initially, we related p a + p ra to p GS by 

log[*(A) - Pr(X) - t Pw ( A)] = log[a(A)] + 6(A)log[/» flJ (A)], (2) 

where a(A) and 6(A) are the function of wavelength and the geometry. [For details, see the MODIS 
water-leaving radiance ATBD Version 2.1 For a given aerosol model, we calculated the a and 6 from 
a linear fit of log(/?t — p r - tp™) to log (p as ) for each wavelength and the given geometry, i.e., for all 
6 vy 00, 4>v and A, from simulations, using the method of least-squares. This resulted in a relationship 
that depends only on the aerosol model and geometry. It reduces the data storage of 3 ~ 4 times 
since the values of p t — p T — tp w and p cs need not be stored for each aerosol concentration, i.e., 
optical thickness (r a ). For the angle 6 V) however, we still have to store a(A) and 6(A) for a large 
number of values of the angle p v for interpolation. To further reduce storage, we expanded a(A) 
and 6(A) in a Fourier series in 6 V and stored only the Fourier coefficients. All of the reflectances 
are even functions of the relative azimuth angle p vy so p t — p r — tp w and p as are as well, and a(A) 
and 6(A) will be even functions of p v . Thus, we can write 

M 

a(0 v ,9, 0,\,4>v) = a (0) (*v,fl 0 ,A)t 2 ^ a< m) (0 v , 9 0 , A) cos m? v , 

rra = 1 

M 

b(9 v ,Bo,\,<f> v ) = bW{0 v ,6 0 ,A)+ 2 b^ m \e v ,9 0 ,X)cos mcf> v , 

m = 1 

0r 

1 f* 

af- m \9 v , 9 0 , A) = - / a(0 v , 0 O , X,<f> v ) cos m<f> v d<t> v , 

* Jo 

1 f T 

b^ m \9 v , 9 0 , A) = - / b(0 v ,0 o ,X,d> v )cos m4> v d$ v . 
r Jo 

a(A) and 6(A) for any <p v can then be obtained easily from the stored values of and 6^ m ^. We 
have now investigated the use of azimuthal interpolation in place of the Fourier transform. This 
will increase the size of the tables; however, it may increase the speed. Also, to see if the accuracy 
of the lookup tables could be improved, we have also replaced Eq. (2) by 

l 0 g[/>t(*) - Pr(A) - *p w (A)] = log [a(A)] + t(^) lo g[Pa>(^)] + c(A)log J (•*)]• (3) 

tMmm mm s 

(if? mm 
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In the implementation of this using the Fourier transform, this increases the size of the lookup 
tables by 50%. 

We have compared the accuracy of the atmospheric correction for standard test cases used in 
the MODIS water-leaving radiance ATBD, i.e., the Maritime, Coastal, and Tropospheric aerosol 
models with a relative humidity of 80%, r a = 0.2 and 0.4, solar zenith angles from 9 0 = 0 to 60°, 
azimuth angle 4> v = 90°, and viewing angle 9 V ^ 45 5 (scan edge). The scan center is not examined 
because the reflectance there is independent of 6 V . The results are presented in Tables 1-6. In 
the tables, the terms “linear” and “quad” refer to Eq. (2) and (3), respectively, the term “Fourier” 
refers to computations with 15 Fourier coefficients for a, b , and c (for the “quad” cases), and the 
numbers, 19, 24, and 37, refer to the number of values of d v used in the azimuthal interpolation. 
Noting that the maximum error in water-leaving reflectance at 443 nm is ±0.001-0. 002, we note 
that there is little difference between the %-arious schemes, with the exception that the quadratic 
fitting procedure [Eq. (3)] produces superior results for the Tropospheric cases with large 9 0 . Thus, 
the choice probably should be based on the speed of the various algorithms. 

(v) To simulate radiative transfer in the 0; absorption, we modified our fifty-layer simulation 
code so that line-by-line computations could be carried out. What we wanted to learn was the 
influence of the O 2 absorption on the radiance or the reflectance leaving the top of the atmosphere 
(TOA) in a band extending from 745 to 7S5 nm (SeaWiFS band 7). Specifically, since the proposed 
SeaWiFS/MODIS algorithm ignores the O 2 absorption, our goal was to be able to estimate w'hat 
the TOA radiance would be in the absence of the absorption band. 

This study has been completed. Two simple modifications to the SeaWiFS algorithm suffice to 
correct for the absorption in the absence of stratospheric aerosols. A possible strategy for detecting 
stratospheric aerosols and/or thin cirrus clouds was also suggested. The complete study is included 
as Appendex 1. A paper describing the results is in press in Applied Optics. 

c. Data/Analysis/Interpretation: See item b above. 


d. Anticipated Future Actions: 
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Table 1: Error in retrieved water leaving reflectance 
for viewing at edge through Tropospheric aerosols. 
True value of r a (S65) is 0.20. 


Algorithm 


0 

0 


0” 

20 ° 

40° 

60° 

linear Fourier 

4-1.2787E-04 

4-] .S096E-04 

+5.7206E-05 

-1.5225E-03 

linear 19 

-h 3.2/8 / E- 0 4 

4-2. 1100E-04 

+ 1.3915E-04 

— 1.4221E-03 

quad. Fourier 

-9.9469E-05 

— 6.1 776E-05 

+ 1.7423E-05 

+ 1.94S1E-04 

quad, linear 19 

— 9.9469E-05 

-3.3934E-05 

+ 1.0664E-04 

+2.8400E-04 

quad, linear 24 

-9.9469E-05 

-6.3206E-04 

+2.8372E-05 

+ 1.8390E-04 

quad, linear 37 

-9.9469E-05 

-5.8398E-04 

+4.6089E-05 

+2.0859E-04 


Table 2: Error in retrieved water leaving reflectance 
for viewing at edge through Coastal aerosols. 
True value of r a (S65) is 0.20. 


Algorithm 


0 

0 


0° 

20° 

40° 

60° 

linear Fourier 

4-2. 1622E-05 

4-1.S054E-04 

4-5.7206E-04 

+ 1.6048E-03 

linear 19 

— 2.1S22E-05 

+ 1.9103E-04 

4-5.9777E-04 

4-1.5984E-03 

quad. Fourier 

-1.4509E-04 

-1.7202E-05 

4-2.9S35E-04 

4-1.2792E-03 

quad, linear 19 

— 1.4609E-04 

— 4.1556E-06 

4-3.2 136E-05 

4-1.2924E-03 

quad, linear 24 

-1.4609E-04 

— 1.7000E-04 

4-2.8661E-04 

4-1.2534E-03 

quad, linear 37 

-1.4609E-04 

' — 1.3556E-05 

+2.9393E-04 

4-1.2632E-03 


Table 3: Error in retrieved water leaving reflectance 
for viewing at edge through Maritime aerosols. 
True value of r a (865) is 0.20. 


Algorithm 

00 

0° 

20 c 

40° 

60° 

linear Fourier 

+6.4034E-05 

4-7.1921E-05 

+5.9443E-05 

— 1.3646E-05 

linear 19 

4-6.3950E-07 

4-7.5364E-05 

+3.0333E-05 

4-7.6612E-05 

quad. Fourier 

— 5.2648E-05 

— 8.71S4E-05 

— 2.1609E-04 

-3.5412E-04 

quad, linear 19 

— 5.2648E-05 

— 7.2438E-05 

-2.0611E-05 

-3.1127E-04 

quad, linear 24 

— 5.2648E-05 

-8.8341E-05 

— 2.1783E-04 

-3.6232E-04 

quad, linear 37 

| -5.2648E-05 

— 8.6554E-05 

-2.1681E-04 

— 3.5245E-04 
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Table 4: Error in retrieved water leaving reflectance 
for viewing at edge through Tropospheric aerosols. 
True value of r a (865) is 0.40. 


Algorithm 


0 

0 


0' 

o 

O 

Cl 

o 

O 

60° 

linear Fourier 

— 1.6432E-03 

-2.069 3 E-03 

— 4.2S73E-03 

— 2. 2461E-02 

linear 19 

— 3.6432E-03 

-2.0169E-03 

— 4.2456E-03 

-1.2328E-02 

quad. Fourier 

- 3 . 246SE-03 

-2.5905E-03 

— 3.4072E-03 

- 6.7699E-03 

quad, linear 19 

— 1.246SE-03 

— 3.54 3 S E-03 

— 3.2235E-03 

-6.6015E-03 

quad, linear 24 

- 2.3468E-03 

— 2.6033E-03 

— 3.3762E-03 

-6.8075E-03 

quad, linear 37 

— 1.1468E-03 

-1.5905E-03 

-3.3365E-03 

-6.7562E-03 


Table 5: Error in retrieved water leaving reflectance 
for viewing at edge through Coastal aerosols. 
True value of r fl (865) is 0.40. 


Algorithm 

0o 

0” 

20° 

40 6 

60° 

linear Fourier 

-2.4262S-04 

+ 2.6253E-04 

H-7.5737E-04 

+ 2.422 3E-03 

linear 19 

-1.4262Z-04 

H- 1 . 8523E-04 

4-8.0943E-04 

+ 1.4326E-03 

quad. Fourier 

-2.1S60E-04 

4-2.9393E-06 

4-4.7064E-04 

+ 2.5163E-03 

quad, linear 19 

— 2.3960E-04 

+2.2 275E-05 

| 

-r5.3098E-04 

“r 2 .oo94E-03 

quad, linear 24 

-2.2960E-04 

4-5.6997E-07 

+4.4502E-04 

+ 2.4777E-03 

quad, linear 37 

-2.2960E-04 ! 

' +6.9279E-06 

-r4.5546E-04 

+ 3.4992E-03 


Table 6: Error in retrieved water leaving reflectance 
for viewing at edge through Maritime aerosols. 
True value of r a (865) is 0.40. 


Algorithm 



0 


0° 

20° 

40“ 

60° 

linear Fourier 

-6.6683E-07 

-4.5924E-05 

-2.6444E-04 

+3.8958E-04 

linear 19 

— 8.4192E-07 

— 4.5866E-05 

— 8.9504E-05 

-4.3686E-04 

quad. Fourier 

-3.7942E-05 

— 9.6642E-05 

— 1.6523E-04 

+ 3.7413E-04 

quad, linear 3 9 

-3.7942E-05 

— 8.4076E-05 

-7.0274E-05 

+5.2993E-04 

quad, linear 24 

-3.7942E-05 

— 2.0493E-04 

-1.2546E-04 

+4.6472E-04 

quad, linear 37 

— 3.7942E-05 

— 2.0312E-04 

-1.2232E-04 

+4.7842E-04 
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(i) We will continue our analysis of the existing simulations from the three-layer code to try 
to understand how to utilize the 1380 nm MODIS band to eliminate the influence of stratospheric 
aerosol. New simulations will be carried out for a variety of stratospheric aerosol types and amounts 
to understand the limits of our conclusions described in Section l.b.(i). 

(ii) We will exercise our fifty-layer MODIS simulation code to provide a better understanding 
of the influence of the aerosol’s vertical profile (in concentration and in type) on the reflectance at 
the top of the atmosphere. In this way we will determine the required inputs into the three-layer 
radiative transfer code for the preparation of atmospheric correction lookup tables. This will start 
when we complete the addition of polarization to the code. 

(iii) We will complete work on the addition of polarization to our MODIS simulation radiative 
transfer code. Our goal is to have a validated code during the first quarter of 1995. 

(iv) We have supplied R. Evans with a new set of lookup tables and access subroutines to see 
the effect of azimuthal interpolation (as opposed to Fourier transformation) on the operation of the 
SeaWiFS — the MODIS prototype — atmospheric correction code. This completes our part of the 
this objective. 

(v) Objective completed. 

e. Problems/Corrective Actions: 

(i) None. 

(ii) Our initial simulations using the MODIS simulation code suggests that a single set of 
lookup tables based on an assumed simple vertical profile for the aerosol may not be sufficient 
to process the data. For example, African dust carried over the Tropical Atlantic is not in the 
boundary layer as assumed in the lookup tables; a separate set of tables must be used to deal w r ith 
such situations. Similar problems arise from the presence of strongly absorbing aerosols, w r hich 
are not in the basic set of models utilized in the present tables, and from the presence of volcanic 
aerosol in the stratosphere, which are absent as well. Thus, the inescapable conclusion is that the 
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lookup tables will have to be regenerated many times using different assumptions in order to be 
prepared to deal with most of the scenarios known to exist. 

Regeneration of the lookup tables is a formidable task. The individual simulations require 
approximately 10 minutes on the DEC 3000/400. Thus, 33,000 simulations would require approxi- 
mately 5500 hours, or approximately 230 days (7.5 months). Utilizing our original four CPU’s on 
a 24-hour basis would reduce the computation time to about 2 months. This was unacceptable for 
many reasons. First if the computers are being used for routine lookup table generation no other 
research can be carried out. Next, the lookup tables will have to be regenerated many times not just 
once, so the computers would rarely be available for other aspects of the algorithm development. 
Finally, after launch of SeaWiFS and processing of the initial data is complete, we will be required 
to generate the necessary new' sets of tables in a very timely manner, i.e., 1-2 weeks not several 
months. We decided we could meet this need with two DECserver 2100 systems, one of w’hich w^as 
procured this year (1994) and with the remaining one scheduled for 1995. 

The DECserver 2100 system has four processors each of w^hich is approximately twice as fast 
as the DEC 3000/400 CPU. Thus, with one DECserver 2100 a single set of lookup tables could 
be generated in one month. With the addition of a second DECserver 2100, the generation time 
w'ould be cut to two weeks. As we plan regular CPU upgrades of the DECservers, by the launch of 
MODIS the speed should be considerably faster which will also allow us to respond to unforseen 
events in a timely manner. Under this plan the four existing DEC 3000/400 w r ork stations will 
always be available for basic algorithm studies, as will the servers, when not generating tables. The 
DECserver 2100 ’s will form the backbone of the Team Member’s Computing Facility to be used for 
post-launch MODIS products. 

(iii) None. 

(iv) None. 

(v) None. 

f. Publications: 
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K. Ding and H.R. Gordon, Atmospheric correction of ocean color sensors: Effects of earth 
curvature, Applied Optics , 33, 7096-7016 (1994). 

K. Ding and H.R. Gordon, Analysis of the influence of 0 2 “A” band absorption on atmospheric 
correction of ocean color imagery, Applied Optics (In press). 

M. Wang and H.R. Gordon, A Simple, Moderately Accurate, Atmospheric Correction Algo- 
rithm for SeaWiFS, Remote Sensing of Environment (In press). 

2. Whitecap Correction Algorithm 

a. Task Objectives: 

As discussed in the previous report, our objective in this period was to build a new radiometer 
to measure the whitecap enhanced reflectance from the sea surface. This radiometer would be 
mounted on a boom off the bow of a ship, directed straight down. The radiometer would record the 
reflectance from a spot of the order of 10 cm in diameter on the sea surface in a continuous manner, 
and a video camera would be bore-sighted with this radiometer to help with data interpretation. 
The background return from clear water would provide the whitecap-free reflectance, and the 
reflectance above the background would provide the average reflectance enhancement of the ocean 
due to whitecaps. A self-contained meteorology package would determine the speed of the ship 
relative to the air, a GPS unit would provide the speed of the ship, and from the two the true wind 
vector would be recorded. Thus the reflectance increase as a function of the wind speed would be 
obtained for algorithm validation. We expected to be able to complete this instrument package by 
the end of the calendar year. 

b. Work Accomplished: 

During the past reporting period we built the whitecap radiometer and associated boom. The 
radiometer measured the upwelling radiance in six spectral bands (as opposed to the proposed 
three). These bands were spread throughout the visible from 410 nm to 860 nm (410 nm, 500 nm, 
550 nm, 670 nm, 760 nm, 860 nm). This system was field tested during a three week MODIS cruise 
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off of Hawaii (R/V Mauna Wave, Chief Scientist Dennis Clark), and the radiometer performed very 
well. To convert the upwelling radiance into reflectance, downwelling irradiance was obtained with 
an old instrument that vre had available (5 channel deck cell). This instrument did not perform as 
well, and needs to be updated. Because of some delays we were not able to obtain the bore-sighted 
video camera or meteorology-station. Fortunately in all these cases redundant data was available 
from other cruise participants. Weather records were recorded by Dr. Stan Hooker and Mr. Jim 
Brown, a video camera was borrowed for short periods from Dennis Clark, and we did side by 
side irradiance measurements with another deck cell operated by Dennis Clark. With all of these 
measurements a very good initial trial data set was obtained. We are currently reducing this data, 
with particular attention to getting correct irradiance measurements for accurate reflectance values. 

c. Data/ Analysis /Interpretation: 

During the cruise mentioned in the previous section a good data set for whitecaps was obtained. 
The wind was typically 20 knots almost, every day, with plenty of whitecaps visible. We are currently 
reducing this data, 

d. Anticipated Future Actions: 

The following improvements or augmentations to the system will be performed in the next 6 
month period. 


1. A video system will be obtained and added to the system to allow 
easier data interpretation. During the sea test of the instrumentation 
it w*as obvious that this is an important addition to the system. 


2. The five channel deck cell will be rebuilt to increase the data to 6 
channels, and to correct some of the stability and reliability problems 
experienced during the field test. 


3. The meteorology package will be integrated into the system. 
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4. The system will be field tested on cruises in which we participate. 

Other actions are to work and complete the data reduction for the initial data set which we 
have already collected. 

e. Problems/Corrective Actions: See item d above. 

f. Publications: 

H.R. Gordon and M. Wang, Influence of Oceanic Whitecaps on Atmospheric Correction of 
SeaWiFS, Applied Optics , 33, 7754-7763 (1994). 

3. In-water Radiance Distribution. 

a. Task Objectives: 

Acquire radiance data at sea. 

b. Work Accomplished: 

Radiance camera is now basically complete. 

c. Data/Analysis/Interpretation: None. 

d. Anticipated Future Actions: 

Acquire data at sea at the earliest opportunity. This will most likely be cruises scheduled by 
Dennis Clark in July 1995 and Fall 1995. 

e. Problems/Corrective Actions: None. 

f. Publications: None. 

4. Residual Instrument Polarization. 
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a. Task Objectives: None. 

5. Direct Sun Glint Correction, 

a. Task Objectives: None. 


6. Prelaunch Atmospheric Correction Validation, 

a. Task Objectives: 

The long-term objectives of this task are two fold. First, we need to demonstrate that our 
atmospheric correction scheme will work to the required accuracy. To effect this we will apply the 
algorithm to compute the sky radiance in the blue from measurements in the near infrared. We 
should be able to do this to about the same accuracy as looking downward from space. Second, 
we need to study the properties of aerosols over the ocean, in particular the aerosol phase function 
and its spectral variation, in order to verify the applicability of the aerosol models on which we are 
basing the atmospheric correction. 

To carry out these objectives requires instrumentation for measuring the sky radiance, solar 
aureole, and the optical thickness of the atmosphere. Instrumentation for measuring the optical 
thickness and the sky radiance is available in our laboratory and has been modified to operate 
with the relevant MODIS spectral bands. Development of an aureole camera for use on a ship 
is a high priority because the sky radiance near the sun is required in aerosol optical property 
retrieval algorithms. Our near-term objective was to see if such inversions could be carried out 
without aureole data. A second objective was to develop and assemble an aureole camera system 
for obtaining such data at sea. The third objective is to begin to obtain a long-term time series 
of the aerosol properties in a maritime environment. For this we are in the process of procuring a 
CIMEL sun/sky radiometer that can be operated in a remote environment and send data back to 
the laboratory via a satellite link. These are similar to the radiometers used by B. Holben and Y. 
Kaufman. 
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b. Work Accomplished: 

We planned and executed a field trip to Key West, FL in April (1994) to obtain sky radiance 
and other aerosol data. Key West is a small island and was chosen to simulate as much as possible 
measurements made from a ship, i.e., to reduce the perturbation by the land as much as possible. 
The planned measurements included (1) aerosol optical thickness as a function of wavelength in 
nine bands from 380 to 1026 nm, (2) sky radiance using the full-sky camera at 560, 671, and 860 
nm, and (3) sky radiance at selected positions using a hand-held radiometer at 558, 669, and 866 
nm. We began analysis of the Key West data; however, for the only day with a clear enough sky, 
the data w’ere obtained at such a small solar zenith angle that inversion is only accurate for the 
forward part of the phase function. Our present inversions yield a value of single scatter albedo 
that was far too low. We are trying to track down the source of the difficulty; however, we believe 
that the error results from the absence of solar aureole data. 

The aureole camera system has now been designed and assembled. The idea is fairly simple, a 
cooled CCD camera has been outfitted with a lens focused at infinity, an interference filter which 
selects the spectral band of interest, and a small aperture over this filter. A neutral density occulter 
is placed approximately 1 meter from the aperture, and this occulter is positioned to block the direct 
sunlight from entering the camera system. In this manner the sky radiance close to the sun can be 
measured, without the direct sunlight (approximately 6 orders of magnitude larger) interfering with 
the measurement. We are currently doing the characterization of the camera system electronics 
and optics. This is part of the calibration procedure. During the next reporting period we will 
also be assembling the field data acquisition software package. Our goal is to field test this system 
during a short cruise in July. 

The CIMEL Electronique, Automatic Sun Tracking Photometer, a Vitel Inc. GOES Data 
Transmitter with accessories, and a set of Solar Power Panels are now on order. 

c. Data/Analysis/Interpretation: Nothing concrete yet. 
d* Anticipated Future Actions: 
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During the next reporting period we will be working on characterizing the camera system, 
and finalizing the data acquisition software. We had also expected to participate on a short cruise 
off of Hawaii during February on which w T e could have tested this system at sea and acquired a 
more complete sky radiance data set. The cruise track for this cruise was optimal for us as it was 
planned to stay in the lee of Lanai, an area where cloud free skies are probable. Due to the slip in 
the SeaWiFS launch, this cruise has now been rescheduled for July. 

We also expect to take delivery of the CIMEL Electronique, Automatic Sun Tracking Pho- 
tometer and assemble to radiometer with the GOES Data Transmitter and Solar Power Panels. 
Testing will begin after assembly with the goal of deploying the instrument in the field in the first 
six months of 1995. 

e. Problems/Corrective Actions: None. 

f. Publications: None. 

7. Detached Coccolith Algorithm and Post Launch Studies, 

a. Task Objectives: 

The algorithm for retrieval of the detached coccolith concentration from the coccolithophorid, 
E. huxleyi is described in detail in our ATBD. The key is quantification of the backscattering 
coefficient of the detached coccoliths. Our earlier studies showed that calcite*specific backscatter 
coefficient was less variable than coccolith-specific backscatter coefficient, and this would be more 
scientifically meaningful for future science that will be performed with this algorithm. The variance 
of the calcite-specific backscatter has been analyzed for only a few species, thus, we need to examine 
this in other laboratory cultures and field samples. There is also a relationship between the rate of 
growth of the calcifying algae and the rate of production and detachment of the coccoliths w r hich 
needs to be further quantified. With this in mind, the objectives of our coccolith studies are, 
under conditions of controlled growth of coccolithophores (using chemostats), to define the effect 
of growth rate on: 
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• the rate that coccoliths detach from cells (which also is a function 
of turbulence and physical shear); 

• the rates of coccolith production; 

• the morphology of coccoliths; and 

• the volume scattering and backscatter of coccoliths. 

The last aspect of these studies will be to perform shipboard measurements of suspended cal- 
cite and estimate its optical backscatter as validation of the laboratory measurements. A thorough 
understanding of these growth-related properties will provide the basis for a generic suspended 
calcite algorithm. As with algorithms for chlorophyll, and primary productivity, the natural vari- 
ance between growth related parameters and optical properties needs to be understood before the 
accuracy of the algorithm can be determined. 

b. Work Accomplished: 

In the last 6 months, we have performed experiments on the effect of turbulence on the rate 
of detachment of coccoliths. These studies were performed under highly controlled conditions of 
turbulent dissipation rates. The rationale for this work was that in order to grow these organisms, 
they must be stirred and since one of our goals is to relate coccolith detachment to growth rate, 
we must insure that we understand the independent effects of turbulence on detachment. 

For the studies of coccolith morphology, we have recently modified our microscope system with 
a phase contrast optical system wdiich gives us much more resolution than we have ever had. This 
will save us considerable time in estimating numbers and morphologies of coccoliths because before, 
we needed to use a scanning electron microscope and make measurements on dead specimens. Now, 
we can count coccoliths and check morphology on living specimens using light microscopy (which 
is essentially instantaneous). In fact, we can actually count the numbers of these 2 micrometer 
diameter plates on individual cells. 
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Successful algorithm development also requires sea-truth data. We continued our analyses 
of optical data from a coccolithophore bloom in the North Atlantic. These data include volume 
scatter measurements, absorption measurements, pigment measurements, cell counts, suspended 
calcite measurements and particulate organic carbon and nitrogen measurements. The first paper 
on this material has been submitted for publication. 

c. Data/Analysis/Interpretation: 

Our results from the laboratory experiments showed that when the cells are actively growing, 
coccolith detachment is not a function of turbulence. However, when cell growth slow T s or stops, then 
plates are held more loosely such that turbulence and coccolith detachment are directly related. 
This represents the first physiological study directly relating turbulence to detachment of coccoliths. 

In regards to our work-up of scatter and absorption data from a coccolithophore bloom in the 
north Atlantic, we have found that chlorophyll-specific absorption by field populations of E. huxleyi 
was similar to previously reported laboratory culture measurements of the same species. Volume 
scatter measurements showed that the suspended coccoliths w T ere responsible for about 80% of the 
total backscatter, 6^ > in the center of the bloom. Vertical profiles of backscatter showed that greatest 
calcite-dependent light scatter was observed just below the base of the mixed layer. Areal maps 
of calcite-dependent backscatter and reflectance (calculated from absorption and backscatter) were 
extremely similar, due to the dominance of backscatter over absorption. Calculated reflectance in 
this mesoscale feature reached as high as 21% at 440 nm and 24% at 550 nm, slightly less than 
what has been observed previously in Gulf of Maine blooms. Total scatter, 6, was also calculated as 
the difference between beam attenuation and absorption. The ratio of backscatter to total scatter 
b h = b b /b was about 0.01-0.02 at 440 nm and 550 nm for b values of 1 to 3 m" 1 . As total scatter 
decreased from 1 bj, increased. The behavior of b b is compared for coccolith- dominated versus 

chlorophyll- dominated waters. Vertical profiles of calcite-dependent scatter, combined with satellite 
remote sensing data, were used to assess the vertical transport velocity of coccoliths (1-2 m per 
day). This velocity was an order of magnitude greater than that expected from theoretical Stokes 
sinking calculations and previous sinking rate experiments. A number of hypotheses are proposed 
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to explain the the discrepancy; the most likely of which is that the subsurface peak in calcite- 
dependent scatter did not result from detached coccoliths sinking but resulted either from plated 
coccolithophores sinking, then detaching their plates, or from deep coccolithophores producing and 
detaching their plates in situ. 

d. Anticipated Future Actions: 

Given the completion of our turbulence experiments, we can now begin the chemostat exper- 
iments in which we will measure the light scatter properties of coccolithophores as a function of 
growth rate. These experiments should start in mid- February. 

There is still more analysis to be done on the field coccolithophore data from the North Atlantic. 
There is sufficient material to publish a second paper on the relationship between backscatter, sus- 
pended calcite concentration, and reflectance. This will be written for Limnology and Oceanography. 

e. Problems/Corrective Actions: None 

f. Publications: None 

8. Other Developments. 

The PI spent a significant portion of the first quarter of this reporting period completing the 
necessary requests to obtain government permission to procure several items of capital equipment. 
Also, all personnel on the project devoted most of their effort in late September and October tow r ard 
revising the Algorithm Theoretical Basis Document (ATBD) for normalized water leaving radiance. 
Version 2 of the ABTD’s for Normalized water-leaving radiance (along with aerosol products) and 
detached coccoliths were delivered to M. King on November 2, 1994. 

The SeaWiFS/MODIS algorithm for estimating the aerosol optical thickness for has now been 
coded and will be delivered to R. Evans in January for incorporation into the processing system. 

A method for combining high- altitude aircraft radiance (upwelling) and surface radiance (down- 
w T elling) for determination of the columnar aerosol optical properties has been developed. A paper 
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on the subject is nearly complete and will be submitted to Applied Optics in January. This could 
provide a powerful method of studying aerosol properties over the ocean. Also, a second study 
concerning the perturbation of the sky radiance measurements made from islands, caused by the 
presence of the island itself, has been carried out and a paper is in preparation. Both of these 
have relevance to the “Prelaunch Atmospheric Correction Validation” (Topic 6 above) portion of 
our research, as well as to the validation of retrieved aerosol properties over the oceans from EOS 
sensors. 
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Abstract 


Two satellite borne ocean color sensors scheduled for launch in the mid 1990’s each have a 
spectral band (nominally 745-785 nm) that completely encompasses the 0 2 A band at 762 nm. 
These spectral bands are to be used in atmospheric correction of the color imagery by assessing the 
aerosol contribution to the total radiance at the sensor. The effect of the 0 2 band on the radiance 
measured at the satellite is studied using a line-by-line backward Monte Carlo radiative transfer 
code. As expected, if the 0 2 absorption is ignored, unacceptably large errors in the atmospheric 
correction result. The effects of the absorption depend on the vertical profile of the aerosol. By 
assuming an aerosol profile — the base profile — we show that it is possible to remove most of 
the 0 2 absorption effects from atmospheric correction in a simple manner. We also investigate the 
sensitivity of the results to the details of the assumed base profile and find that, with the exception 
of situations in which there are significant quantities of aerosol in the stratosphere, e.g., following 
volcanic eruptions or in the presence of thin cirrus clouds, the quality of the atmospheric correction 
depends only weakly on the base profile. Situations with high concentrations of stratospheric 
aerosol require additional information regarding vertical structure to utilize this spectral band in 
atmospheric correction; however, it should be possible to infer the presence of such aerosol by a 
failure of the atmospheric correction to produce acceptable water -leaving radiance in the red. An 
important feature of our method for removal of the 0 2 absorption effects is that it allows the use of 
lookup tables that can be prepared in the absence of 0 2 absorption using more efficient radiative 
transfer codes. 
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Introduction 


The radiance L w exiting the ocean in the visible carries information concerning the concentra- 
tion of phytoplankton in near-surface waters. Phytoplankton are microscopic plants that through 
photosynthesis, i.e., they combine C0 2 (dissolved in the water) and water to produce carbohy- 
drates, form the first link in the marine food chain. Their influence on the C0 2 makes knowledge 
of their spatial and temporal variability important in the global carbon cycle. The Coastal Zone 
Color Scanner 1 (CZCS) launched by NASA on Nimbus 7 in the fall of 1978 (and operational until 
mid 1986) acquired imagery over the oceans in four spectral bands centered at 443, 520, 550, and 
670 nm, with spectral widths of approximately 20 nm. Gordon et al. 2 showed that the phytoplank- 
ton pigment concentration, C — the sum of the concentrations of chlorophyll a and its degradation 
product phaeophyton a — could be derived from this imagery after removing the effects of the 
atmosphere. The CZCS spectral bands were located in atmospheric "windows,” so the atmospheric 
effects were due almost entirely to scattering by molecules (Rayleigh scattering) and by aerosols. 
Except for a small variation due to changes in atmospheric surface pressure, the Rayleigh scatter- 
ing contribution can be computed exactly. 2 In contrast, the aerosol scattering is highly variable in 
space and time and must be determined from the imagery itself. This was effected by utilizing the 
band at 670 nm, where L w is very small so the top of the atmosphere (TOA) radiance is due to the 
atmosphere, to assess the aerosol’s contribution there. The spectral variation was then estimated 
by either using regions in the imagery where C < 0.25 mg/m 3 (clear water regions 4 for which L w 
is known in the green bands as well as the red band) or by assigning a nominal spectral variation 1 
based on experience and on the properties of maritime aerosols. However, each procedure required 
assumptions, preventing a totally deterministic correction. 

With the success of the CZCS as a proof-of- concept mission several ocean color instruments are 
being prepared for launch: the sea-viewing wide-field-of-view sensor (Sea^ViFS); the moderate res- 
olution imaging spectroradiometer (MODIS); 8 and the ocean color and temperature sensor (OCTS) 
on the advanced earth observation satellite (ADEOS) to be launched by Japan. To provide better 
atmospheric correction, these instruments have spectral bands in the near infrared NIR, where the 
ocean can be considered to be a Fresnel-reflecting medium that absorbes all photons penetrating 
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the surface. In the case of SeaWiFS and OCTS, the NIR bands are 40 nm wide and positioned at 
765 and 865 nm in windows free of water vapor absorption. The rather large spectral widths were 
governed by the requirement for high signal-to-noise ratios. Thus, the 765 nm band on both sensors 
was forced to completely encompass the O 2 "A” absorption band at 762 nm. It was believed that 
correction for the ~ 10% absorption due to the 0 2 “A” band would be possible. We have developed 
an approach to the correction of the 765 nm bands for both SeaWiFS and OCTS for O 2 absorption 
and report such in this paper. 

We begin with a brief review of the proposed 9 atmospheric correction algorithm for SeaWiFS 
which utilizes the 765 nm band, and which ignores the presence of the O 2 absorption, i.e., was 
developed assuming the O 2 absorption feature in the 765 nm band did not exist. Next, through 
a series of Monte Carlo simulations, we examine the influence of the O 2 absorption on radiative 
transfer in the atmosphere. Based on these simulations a method for assessing and removing the 
O 2 effect is developed. Finally, we present examples of the performance of the end-to-end SeaWiFS 
atmospheric correction algorithm in the presence of the O 2 absorption. 

The proposed SeaWiFS correction algorithm 

In a recent paper, Gordon and Wang 9 proposed an algorithm for the atmospheric correction 
of SeaWiFS, and ultimately, MODIS. Briefly, the total radiance L t { A) measured at the top of the 
atmosphere at a wavelength A can be decomposed as follows: 

L t { A) = Lr{ A) + L a { A) + I ro (A) + t{9 v ,\)L w (X), (1) 

where £ r (A) is the radiance resulting from multiple scattering by air molecules (Rayleigh scattering) 
in the absence of aerosols, £ a (A) is the radiance resulting from multiple scattering by aerosols in the 
absence of the air, and £ ra (A) is the interaction term between molecular and aerosol scattering. 10 In 
this equation, t is the diffuse transmittance of the atmosphere along the viewing direction specified 
by 9 V , the angle between the normal to the sea surface and the sensor. 11 Radiance arising from 
specular reflection of direct sunlight from the sea surface and subsequent propagation to the top of 
the atmsophere without scattering (sun glitter) has been ignored because SeaWiFS can be tilted 
away from the glitter pattern. However, specular reflection of skylight, created by both molecular 
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and aerosol scattering, is included in L r , L ai and L ra , The influence of whitecaps has been ignored 
under the assumption that their contribution can be estimated from an estimate of the surface 
wind speed. 12 

The goal of the atmospheric correction is the retrieval of L w from L t . It is convenient to 
convert radiance (L) to a reflectance ( p ) defined to be 7ri//o cos #o, where F 0 is the extraterrestrial 
solar irradiance, and 8q is the solar zenith angle. With this definition, Eq. (1) becomes 

p t (A) = p r (A) + Pa(A) + p ro (A) + t(^ v > (2) 

The reason for the particular decomposition in Eqs. (1) and (2) is that p r can be computed exactly 
given the surface atmospheric pressure and A. 3 

After computation of p r , the algorithm utilizes p t — p r = p a 4* p T a in the near infrared (NEEt) 
at 765 and 865 nm, where p w can be taken to be zero except in turbid coastal waters or possibly 
intense open-ocean coccolithophore blooms, 13 to choose two aerosol models from a set of candidate 
models. The candidate models were taken from those proposed by Shettle and Fenn. 14 The chosen 
aerosol models are then used to predict p Q + p ra in the visible which, when combined w r ith p r 
there, yields the desired tp w . To assess the efficacy of the algorithm, simulations of p t using aerosol 
models which were similar, but not identical, to the candidate aerosols models, were carried out. 
The simulated p t was then used as pseudo data for insertion into the correction algorithm. The 
error A p in the retrieval of tp w at 443 nm was nearly always found to be < 0.002 and often 
< 0.001. This error meets the SeaWiFS goal of retrieving L w (or p w ) at 443 nm to within ** ±5% 
in waters wfith low phytoplankton pigment concentrations, e.g., the Sargasso Sea in summer. 

The algorithm as described ignores the presence of the 02, he., in the simulations of p t and 
p r the model atmospheres were free of the O 2 absorption in the 765 nm band. In the presence of 
O 2 absorption p t will be smaller in the 765 nm band than in its absence, so the computed value of 
Pa + Pra &t 765 nm will be smaller. Operation of the present correction algorithm in the presence 
of the O 2 absorption with the smaller value of p a + p ra will lead to the choice of an inappropriate 
model with which to determine p a + p ra in the visible. 
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Our approach to dealing with the effects of the O 2 absorption in the 765 nm band is to use the 
value of p a + p ra derived in the presence of the O 2 absorption to determine the value it would have 
were the O 2 absorption absent. The simulations presented below suggest that such an approach is 
viable. 

Radiative transfer In the 0 2 A band 

In this section we describe the technique used to study the radiative transfer in the O 2 A 
band. What we wish to learn is the influence of the O 2 absorption on the reflectance p t leaving the 
TOA in a band extending from 745 to 785 nm. Specifically, since the proposed SeaWiFS algorithm 
ignores the O 2 absorption, our goal is to be able to estimate what p a -f- p ra would be in the absence 
of the absorption band. 

A. The O 2 A band. 

The O 2 A band extends from about 759 to 770 nm. There are 286 individual absorption lines 
having appreciable line strengths. Figure 1 shows a sample of the sea-level absorption coefficient as 
a function of wave number, u = 1/A, computed from the positions, line strengths, and line widths 
of the individual absorption lines (assuming a Lorentz shape) taken from the AFGL compilation. 15 
Note that this figure covers only about 1.165 nm but contains 12 spectral lines, and the absorption 
coefficient varies from approximately 30 km^ 1 to about 0.2 km** 1 over a small fraction of a nm. 
Because of the strong variation of the absorption coefficient with wavelength, the absorption over a 
band containing several lines wull not be an exponential function of the path length. Thus, it is not 
possible to assign a single mean absorption coefficient to the entire O 2 A band. Furthermore, since 
the individual spectral lines are pressure and temperature broadened, even at discrete wavelengths, 
i.e., bands with wudth <C the wddth of the individual spectral lines, the absorption coefficient will 
be dependent on altitude in the atmosphere. Thus, a complete treatment of the radiative transfer 
in this absorption band requires an atmosphere consisting of several layers in w T hich the absorption 
coefficient is a very strong function of frequency. 
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B. Preliminary discussion of the 0 2 effect. 

It is possible to understand qualitatively the effect of the 0 2 A band on the radiance exiting 
the top of the atmosphere by examining single scattering. First, we assume that the atmosphere is 
free of aerosols, i.e., we only have Rayleigh scattering and the scattering coefficient will vary with 
altitude in proportion to the density in the same manner as the 0 2 abundance. For a given viewing 
geometry we define the air mass M as 

M = — H 

COS V v COS Uq 

where 9 0 is the solar zenith angle and 9 V is the viewing angle, i.e., the angle between the surface 
normal and the direction of propagation of the radiance exiting the TOA. Photons scattering from 
molecules at any altitude will have traversed a path of length proportional to M upon exiting the 
atmosphere. Thus, w r e expect the decrease in radiance exiting the atmosphere to be a function of 
M\ albeit not an exponential function. In the case of multiple scattering the path of the photon 
is no longer proportional to M so a similar argument does not apply; however, since the Rayleigh 
scattering optical thickness is small 0.025) at 765 nm, multiple Rayleigh scattering w ? ill be small 
and the radiance decrease will still depend on M in much the same manner as for single scattering. 

The addition of aerosols causes tw^o complications: the aerosol concentration is a strong and 
variable function of altitude; and the aerosol concentration is usually sufficiently high that multiple 
scattering is significant. The influence of the vertical profile of the aerosol concentration is easy to 
understand in the single scattering approximation. Typically, over the oceans most of the aerosol is 
in the marine boundary layer which is 1-2 km thick. The aerosol component of the TOA radiance 
at 765 nm with a high concentration of aerosol in the boundary layer will be significantly larger 
than the molecular-scattering component. This radiance will have had to travel through most 
of the atmosphere (twice) before reaching the TOA. The total path is proportional to Af, so we 
expect that the radiance decrease will be larger than that for an aerosol-free atmosphere (because 
more photons backscattered to the TOA will have had to travel farther in the atmosphere), and 
that it will be a function of the airmass. In contrast, if there is a high concentration of aerosol 
in the stratosphere, e.g., following a major volcanic eruption, a fraction of the TOA radiance will 
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have scattered from the stratosphere and not have traveled through a significant portion of the 
atmosphere. Then, for the same aerosol concentration as in the marine boundary layer example, 
the fractional decrease in the radiance due to the O 2 absorption will be less. 

In general the aerosol will be distributed continuously throughout the atmosphere. A typical 
profile will display a high concentration in the boundary layer, a rapid decrease in the free tropo- 
sphere and slower decrease in the background stratosphere, and occasionally an increase (sometimes 
large) due to recent (e.g., the past 1-3 years) volcanic activity. 16 Thus, to simulate the fractional 
decrease in the TOA radiance near 765 nm due to the O 2 A band absorption requires an atmosphere 
composed of enough layers to adequately represent the vertical distribution of the aerosol. 

C. The radiative transfer code and model. 

The radiative transfer code we developed for this study is a backward Monte Carlo for a plane 
parallel atmosphere. Polarization of the light field is ignored. The backward procedure is used 
because it provides a significant increase in accuracy (over the forward Monte Carlo), for a given 
expense in computational time, when the radiance is desired in only a single viewing direction. In 
this kind of simulation, photons are ejected from the detector and followed to the sun. Actually, at 
each interaction in the atmosphere the probability that the photon will be scattered to the TOA in 
the direction of the sun is computed and used as the estimator. The atmosphere is bounded by the 
sea, which we model as a Fresnel reflecting surface that absorbs all photons that are transmitted 
through it. This is realistic because of the large absorption coefficient of liquid water in the NIB,. 
The possibility that at any interaction a photon could be scattered toward the sea surface and be 
Fresnel-reflected back to the TOA is included in the Monte Carlo estimator. The code includes 
provision for the sea surface to be roughened by the wind wdth surface slopes obeying the Cox and 
Munk 17 distribution; however, in all of the simulations presented here the surface was smooth. 

For our model atmosphere we assumed that the density and temperature follow the U.S. 1976 
Standard Atmosphere. 18 O 2 is distributed in proportion to the atmospheric density. The individual 
O 2 spectral lines are assumed to have a Lorentz line shape. The line strengths and line widths 
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(at NTP) were taken from the AFGL 15 compilation, and the line width varies with altitude in 
proportion to P/y/T , where P is the pressure and T is the absolute temperature. The full width at 
half maximum (FWHM) at the sea surface of a typical absorption line in the A band is of the order 
of 0.1 cm“ l or ~ 0.006 nm (Figure 1). The atmosphere is composed of thirty-seven homogeneous 
layers. For each layer the O 2 .absorption coefficient as a function of wave number v was computed 
and tabulated at 10 5 values of v between v — 13,422.82 cm“ l (A = 745 nm) and 12,738.85 cm~ l 
(A = 785 nm). With this resolution, there are approximately 15 tabulated values wuthin the FWHM 
of each absorption line at the surface. Linear interpolation is used for values of v between tabulated 
values. 

For a realistic treatment of the aerosol in the code, we further divide the atmosphere into 
four broad regions (Figure 2): (1) the marine boundary layer from the surface to 2 km, where 
the aerosol concentration is independent of altitude; (2) the free troposphere, where the aerosol 
concentrations varies in proportion to exp [— 2 // 1 ], where 5 is the altitude (2-12 km) and h (called the 
scale height) is 2 km; (3) the background stratosphere (12-30 km), where the aerosol concentration 
is also exponential with a scale height of 5 km; and (4) a volcanic region (18-20 km) within the 
stratosphere that can contain a uniformly mixed volcanic aerosol. The optical properties of each 
of the four regions can be characterized by individual aerosol models, and any of the regions can 
be free of aerosols if desired. There is no aerosol above 30 km. 

As in Ref. 9, the aerosol models used for the two lower regions were taken from Shettle and 
Fenn. 14 Based on size distribution and composition measurements, they developed two models 
called the Tropospheric (to model the aerosol in the free troposphere — few large particles) and 
Oceanic (to model the aerosol produced by sea spray — few small particles). They combined these 
to form the Maritime model to represent the aerosol in the marine boundary layer. In terms of 
total aerosol number per unit volume, the Maritime model consists of 99% Tropospheric and 1% 
Oceanic. Gordon and Wang 9 added a Coastal model (99.5% Tropospheric and 0.5% Oceanic) to 
provide a description of the aerosol that might be more representative of the boundary layer near 
the coast. 
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For the upper two regions we use a model for the background stratosphere from the WMO 
and a model for the volcanic aerosol from King et al. :o Both models assume a 75% solution of 
H2SO4. We have also included a volcanic ash model 19 to represent fresh volcanic aerosol. The size 
distributions in these regions are modified gamma distributions. 

With all of the optical properties of the atmosphere prescribed, the Monte Carlo code operates 
as follows. A photon is ejected from the receiver in a direction opposite to the observation direction 
with a w’ave number v chosen from a uniform probability density over 12,738.82 < v < 13,422.82 
cm -1 . This choice of v corresponds to assuming the extraterrestrial solar irradiance is constant 
over the this spectral interval, i.e., F 0 is replaced by its spectral band-averaged value. This is a 
convenience rather than a necessity. The Rayleigh optical extinction coefficient is then computed 
for v based on Fenn et al. 16 The aerosol extinction coefficient for each layer is provided by the 
aerosol model and concentration. It is assumed to be independent of A over the 745-785 nm 
band. The optical distance (defined as / c(l)dl, where c{l) is the extinction coefficient along the 
path /) the photon travels is then determined from an exponential probability density function. 
Prom the tabulated values of the absorption and scattering coefficients, the physical distance that 
this corresponds to is determined, yielding the position of the first interaction point. The various 
estimators are computed, the photon allowed to scatter, and the process is repeated. In a typical 
computer run we trace the histories of 10^ photons. This implies that there are approximately 
1,300 photons sampling the FWHM of each individual Oj absorption line at the sea surface. 


D. Tests of the radiative transfer code. 

Validation of the radiative transfer code was effected in two ways. In the absence of the Oj 
absorption, the code was compared to the output of a successive order of scattering code similar 
to that used by Gordon and Wang. 9 In this case, v was not sampled, the molecular scattering 
was distributed uniformly from the surface to 100 km (r r = 0.0255), and the aerosol scattering 
was uniformly distributed from the surface to 30 km (r„ = 0.2). The difference in the radiances 
computed by the two codes was usually < 0.1%. This suggests that when the 0 2 absorption 
features of the code are not employed, it is capable of returning radiances accurate to within 0.1%. 
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It is difficult to test completely the code’s handling of the 0 2 absorption band, because, in 
contrast to the situation in the previous paragraph, we do not have a second, highly accurate code 
for comparison. However, the LOWTRAN 21 atmospheric propagation code, which treats band 
absorption well, does have an approximate multiple scattering capability for calculating radiance 
reflected from the earth-atmosphere system. 22 If aerosol scattering is omitted, at 765 nm multiple 
scattering effects should be small enough so that LOWTRAN can compute the radiance. Since 
LOWTRAN does not have the provision for a Fresnel-reflecting sea surface, we ran both codes 
with a totally absorbing lower boundary. The differences in the radiances (i r ) over the spectral 
band 12,735 <v< 13,425 cm“ ] between LOWTRAN and our code were 1%. Since our code 
treats both multiple scattering and gas absorption in a more precise manner than LOWTRAN, we 
believe this demonstrates that our code is satisfactory. 


Finally, since the Lorentz line shape underestimates the absorption in the core of the lines and 
overestimates the absorption in the wings in the upper atmosphere (z > 10 km), we carried out a 
series of simulations using the Voigt 23 Line shape to describe the 0 2 absorption. For an aerosol-free 
atmosphere, the two line shapes yielded reflectances within the inherent error of the Monte Carlo 
procedure 0.1%). The Lorentz line shape was used in the computations reported here. 

Removal of 0 2 A band absorption 


In this section we develop the necessary relationships for removing the effects of the 0 2 A 
band absorption from the 765 nm band. We shall use the notation that primes represent quantities 
computed or measured when the 0 2 absorption band is present ( prime and present both begin with 
“p”), and unprimed quantities represent those in the absence of the 0 2 absorption. Thus Eq. (2) 
becomes 


( 3 ) 


Pt “ Pr — Pa + Pra — PAi 
Pt Pr Pa d* Pra — P A' 

where we have explicitly used the fact that p w is taken to be zero. The known quantities in these 
equations are p[ (measured) and p r . The SeaWiFS correction algorithm 9 will employ precomputed 
lookup tables of p r at standard atmospheric pressure, and p\ for a variety of aerosol models and 
optical thicknesses. As the notation suggests, these tables for the 765 nm SeaWiFS band have been 
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computed assuming that the O 2 absorption band does not exist. Thus, we need to find p\ from p l t 
and p T . 


Our technique involves three steps: (1) finding p' r from an empirical relationship between p' r 
and p r developed from simulations; (2) using the second equation in Eq. (3) to deduce p^\ and 
(3) finding px from an empirical relationship between p\ and p x developed from simulations. The 
desired px is then used in the atmospheric correction algorithm. A desirable feature of this 0 2 
correction is that both sets of 0 2 -free lookup tables can still be employed. This is important 
because of the intense computational burden required for computing accurate reflectances in the 
presence of the 0 2 absorption. In this section the empirical relationships we have mentioned above 
are developed. 

A. Relationship between p T and p' T . 

As described in an earlier section, we expect that the fractional decrease in the reflectance of 
the atmosphere-ocean system in the 745-785 nm spectral band due to the 0 2 absorption will be a 
function of the two-way airmass M . Thus, we carried out a series of Monte Carlo simulations for 
an aerosol- free atmosphere using a variety of solar zenith and viewing angles. Both p T and p' T are 
computed, from which we found the fractional change in reflectance, ( p T — p' r )j p f ri as a function of 
M . The results are provided in Figure 3. Note that the minimum air mass is 2, and for this the 
0 2 absorption results in a decrease in the reflected radiance by approximately 7%. Typical viewing 
geometries have 2 < M < 5, so the decrease in reflectance will vary between - 7 and ~ 11%. The 
dotted line in Figure 3 is a least-squares fit of the computations to the air mass, i.e., 

log 10 (j~ — 7 "^"^ — GrO + a r\M + a r2 M" = P r(M), 

where the a r0 , a rl> and a r2 are -1.3491, +0.1155, and -7.0218 x 10“ 3 , respectively. Using this fit, 
the desired p ! r for any M can be estimated from 


= (1 + io p -<">) "V 


( 4 ) 


Figure 4 provides the error (%) in the estimated value of p f T , using Eq. (4), for the simulations in 
Figure 3. It suggests that p ' T , the reflectance when the 0 2 absorption is present, can be estimated 
from p r , the reflectance in the absence of the 0 2 absorption, to within ^ ±0.2% solely from the 
sun- viewing geometry. 


B. Relationship between p A and p\. 

Now we develop a relationship between p A and We proceed in the same manner as with 
the Rayleigh component, i.e., we carry out a series of Monte Carlo simulations for an atmosphere 
including both aerosol and molecular scattering; however, there is an inherent difficulty in that the 
relationship cannot be independent of the vertical structure of the aerosol. Therefore, the estimate 
of px from p A will be less accurate than the estimat ion of p' r from p r . Clearly, w ? e must choose an 
aerosol profile, on which to base the relationship, that is as representative as possible. We call this 
the “base” profile. 

To assess the magnitude of the error in the relationship w'hen an unrepresentative base profile 
is chosen, we begin w r ith a very simple base profile and find the error in the prediction of p A from 
p l A . For this purpose, the base profile will consist of a uniformly mixed aerosol in the boundary 
layer (thickness of 2 km), with the free troposphere and stratosphere being free of aerosols. Thus 
all of the aerosol is in the marine boundary layer. The base aerosol optical thickness r 0 is taken to 
be 0.2693 and the aerosol model is the Shettle and Fenn 14 Maritime model at 50% relative humidity 
(M50). Figure 5 provides [p A - p' A )/ p A as a function of M for a series of simulations using this 
base aerosol. As in the case of a pure Rayleigh scattering atmosphere, the line in Figure 5 is a 
least-squares fit of the computations to the air mass, i.e., 

log 10 -) = + a Al M + a A2 M 2 = P A (M ), (5) 

where the a’s are constants. Using this fit, the desired p A for any M can be estimated from 

PA= (l + 10^ (M) )p' A . (6) 
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We carried out simulations of p\ for a set of aerosol profiles, concentrations, and type to test the 
accuracy with which p A could be estimated given p[ y if the three steps described at the beginning 
of this section were used. First we estimated p' r from p r using Eq. (4). This provided p f A = p[ — p f T . 
Then we estimated p A from p\ using Eq. (6), based on their relationship in Figure 5, for the simple 
base profile. Figures 6 and 7-report the error, Ap Ai in this prediction of p A from p[. In Figure 7 
the Shettle and Fenn 14 Tropospheric aerosol at 50% relative humidity (T50) has been used as the 
aerosol type in all layers, although the base aerosol was still M50. Two important conclusions can 
be drawn from the computations presented in these figures: (1) when the base profile is correct, 
the error in p A is < 0.5% for M < 5, and nearly independent of the aerosol concentration and 
the aerosol type; and (2) even with the incorrect base profile, the error is usually < 1% as long 
as the test profile is not “too different” from the base profile (the case with a 5 km scale height in 
the free troposphere would be an example of a “too different' 1 profile). This suggests that using the 
methodology described here, p A can be derived from p\ with an error < 1% as long as a realistic 
base aerosol profile is employed. 

C. Choice of the base aerosol profile and type. 

FVom the results presented above it is clear that it is important to use as realistic a base aerosol 
profile as possible. In fact, the computations suggest that the base profile should have some aerosol 
in the free troposphere and the stratosphere (unless these happen to be aerosol free) in order to be 
assured that the error in p A < 1%. With this in mind, we choose the base aerosol in the following 
manner. Referring to Figure 2, we use the Shettle and Fenn 14 Maritime model with 80% relative 
humidity (M80) in the marine boundary layer, and their Tropospheric mocel with relative humidity 
of 50% (T50) for the free troposphere. The background stratospheric model 19 (B) is proposed for 
the stratosphere. The volcanic component is taken to be absent in the base profile. The aerosol 
concentrations are adjusted so that when the visible range (defined 14 to be 3.96 divided by the 
total extinction coefficient at 550 nm) is 25 km, the aerosol extinction coefficient is a continuous 
function of altitude, i.e., there is no discontinuous jump in the aerosol extinction from the top of the 
boundary layer to the bottom of the free troposphere of the type shown on the figure. Note, for any 
wavelength other than 550 nm, discontinuities in aerosol extinction will occur at the boundaries of 


14 


PASS & 

a* mmmmm 


each region because the extinction coefficient for each region will have a different spectral variation. 
The base aerosol is denoted by M80-T50-B, which stands for M80 in the boundary layer, T50 in 
the free troposphere, and B in the stratosphere. Similar notation will be used for aerosol profiles 
through out the remainder of this paper. The (p A - p f A )/ p A - M relationship for this base aerosol 
profile is provided in Figure 8. Comparison with Figure 5 shows that for a given M , (p A - p'a)! p' A 
is slightly smaller (less O 2 absorption) than it is for the base profile with all of the aerosol in the 
marine boundary layer. For the M-80-T50-B base profile, the values of a A\ > ^d a^2 in Eq. (5) 
are -1.0796, +9.0481 X 1CT 2 , and -6.8452 x 10~ 3 , respectively. 

Application to atmospheric correction 

In this section we apply the O 2 absorption correction scheme developed above to simulated 
ocean color data, e.g., SeaWiFS or OCTS. We will examine the optimum case, i.e., assume that 
the aerosol profile in the free troposphere and stratosphere is identical to the base profile. Only 
the aerosol in the boundary layer is allowed to be different from the base profile in concentration 
and/or in type. This of course would require knowing the free troposphere and stratospheric aerosol 
profile in advance. 

We examined the following cases: (1) a visible range ( VR ) of 40 km, 'with the Maritime (M80) 
or Tropospheric (T80) aerosol at 80% relative humidity in the marine boundary layer; and (2) the 
same aerosol models with VR = 15 km. These correspond to either ~ 40% less (case 1) or ^ 72% 
more (case 2) aerosol in the boundary layer than in the base profile. The total and boundary-layer 
aerosol optical thicknesses for these models are provided in Tables 1-4. 

Simulations were carried out for both p[ and p t so we could compare the performance of 
the algorithm in the presence and absence of the O 2 absorption to test the efficacy of the O 2 
correction. The geometry of the simulations consisted of solar zenith angles 9 0 = 0(20°)60°, witb 
Q v = 0° and ^ 45° corresponding, respectively, to viewing at the scan center and the scan edge. All 
the simulations were carried out using the Monte Carlo code described above. The simulated value 
p\ (in the presence of the O 2 absorption) at 765 run was then combined with p T through Eqs. (4) 
and (6) to estimate p A , w r hich was then used in the SeaWiFS correction algorithm. For comparison, 
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in the absence of the O 2 absorption, the simulated pt was used in the normal way in the correction 
algorithm. 9 The measure of effectiveness of the combined algorithm (O 2 removal and atmospheric 
correction) is the ability to recover 443 nm with an error of < 0.002. In this paper, the 

particular implementation of the algorithm provided in Ref. 9 utilizes the Maritime, Coastal, and 
Tropospheric models at relative humidities of 50%, 70%, 90%, and 99% as twelve candidate aerosol 
models. As in Ref. 9, it is assumed in this implementation that the aerosol is confined to the 
marine boundary layer. 

Figures 9a-9d and lOa-lOd show the error in the retrieved value of tp w (called A p on the figures) 
as a function of 6 0 and 9 V , for the four test cases: MS0-T50-B and T80-T50-B with VR = 15 and 
40 km. The open circles represent the error in the retrieval if the spectral band from 745 to 785 nm 
were free of 0 2 absorption. Note that even if the 0 : absorption were absent, A p can at times be 
quite large. This owes to the fact that in the preparation of the lockup tables for implementation 
the SeaWiFS algorithm it was assumed that the aerosol was all in the marine boundary layer, 
i.e., a two-layer radiative transfer code (aerosols on the bottom) was employed. This assumption 
leads to errors when the aerosol is not confined near the surface. Even so, Figure 9 shows that 
when the aerosol optical depth at 443 is < 0.4 the correction error is within the desired range 
range, i.e., ±0.002. In contrast, when the optical thickness at 443 reaches 0. 7-0.8 (Figure 10), the 
error can become excessive. Measurements suggest that optical thicknesses in the latter range are 
not common when the aerosol is locally generated over the oceans.'^’* These simulations suggest 
that, even in the absence of the O 2 absorption, e.g., the atmospheric correction of MODIS, the 
assumptions made by Gordon and Wang 9 may need to be modified by the addition of a third 
intermediate layer containing a mixture of aerosol and Rayleigh scatterers. 

However, our main concern here is the efficacy of the removal of the effects of the O 2 absorption 
from the atmospheric correction. For this purpose we need only compare the differences between 
the open circles and the filled circles in Figures 9 and 10. We note that the difference is usually 
< 0.001 and often considerably less, indicating that the 0 2 correction works reasonably well. 
Recall that if p A is too small, an undercorrection will occur, i.e., the water-leaving reflectance at 
443 will be too large (A p > 0). Thus, in Figure 9, p A is a little too large, while in Figure 10 it is 


too small. This is consistent with the fact that the (p A - p' A )h' A ~ M relationships in Figure 8 
will move upward (more absorption) if the aerosol concentration in the boundary layer is increased 
and downward if the concentration is decreased with respect to the base profile. This implies that 
for the cases in Figure 10, P A {M) in Eq. (6) will lead to a value of p A that is too small, while for 
the cases in Figure 9, the reverse is true. 

To provide an appreciation of the magnitude of the O 2 effect that has been removed, in Table 
5 we present the average A p over 0q for each model when the existence of the O 2 absorption is 
simply ignored and the proposed SeaWiFS algorithm operated using p[. Comparison with the 
results in Figures 9 and 10 suggest that we have removed about 9f% of the error in the worst case 
(T80-T50-B at the scan edge). 

In the simulated test of the O 2 absorption correction above, we used a base profile that was 
correct for the free troposphere and the stratosphere. This should be the optimum situation. To 
try to understand the effect of choosing an incorrect base profile on the end-to-end atmospheric 
correction algorithm, we computed A p using the base profile that was used in the preparation of 
Figure 5, i.e., a base profile assuming that all of the aerosol is in the boundary layer and modeled 
as M50. This is in contrast to the base profile used in the preparation of Figures 9 and 10 in which 
~ 30% of the aerosol optical thickness at 765 nm was above the boundary layer. The resulting 
corrections at the scan center using M50 as the base profile are provided in Figure 11. Comparing 
Figures 9a with 11a, 9b with lib, 10a with 11c, and 10b with lid, we see that the filled circles are 
< 0.001 (and typically much less) lower using the M50 base profile compared to the M80-T50-B 
base profile. Thus, a large error in the base profile does not significantly influence the overall 
atmospheric correction. In fact, for the VR = 15 km case, the incorrect base profile actually 
produces better corrections. Similar results Eire obtained at the scan edge (not shown). 

Volcanic Aerosols 

A situation in which the base profile can strongly influence the atmospheric correction is that 
following volcanic eruptions, e.g., El Chichon, which can inject significant quantities of aerosol into 
the stratosphere. To understand the effect of this on the 0 2 correction presented above, we have 
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examined four cases with volcanic aerosol. We used the El Chichon aerosol model developed by 
King et al. 20 and considered the aerosol to be concentrated in a layer from 18 to 20 km (Figure 
2). The optical thickness at 550 nm was taken to be 0.05 or 0.10. Figure 12 provides the results of 
simulated retrieval errors Ap at the scan center as a function of 0q. The notation for the aerosol 
profile now includes “V” followed by a number indicating the volcanic aerosol optical thickness at 
550 nm, e.g., M80-T50-B-V05 indicates the addition of a volcanic aerosol with optical thickness 0.05 
at 550 nm, etc. Examination of Figure 12 shows that the error in the presence of the 0 2 absorption 
(filled circles) is always below (more negative) the error what would occur in the absence of 0 2 
absorption (open circles). This occurs because the presence of the volcanic aerosol causes less 0 2 
absorption than the base profile. Thus, Pa(M') in Eq. (6) is too large so px w’ill be too large 
resulting in an overcorrection (A p < 0). The effect is more significant the greater the contribution 
from the volcanic aerosol compared to the boundary-layer aerosol. This can be seen by comparing 
Figures 9a, 12a, and 12b for which the volcanic optical thickness progresses from 0 to 0.05 to 0.10 
with V R — 40 km, and Figures 10a, 12c, and 1 2 d , for V R = 15 km. For V R — 40 km, we see a slow 
degradation in the quality of the atmospheric correction even in the absence of 0 2 absorption with 
increasing volcanic optical thickness and a large overcorrection in the presence of 0 2 absorption. 
In contrast, for V R — 15 km we see the same slow degradation in the absence of 0 2 absorption, 
but a smaller 0 2 absorption overcorrection. 

It is dear that application of the proposed 0 2 absorption correction in the presence of volcanic 
aerosol will require information concerning the vertical distribution, concentration, and optical 
properties of the volcanic aerosol. 

Summary and Conclusions 

By simulating the radiance reflected from the ocean-atmosphere system in the 0 2 A absorp- 
tion band using a line-by-line backward Monte Carlo radiative transfer code, we have shown the 
influence of the absorption band on the atmospheric correction of ocean color sensors utilizing a 
~ 745-785 spectral band, e.g., SeaWiFS and OCTS. If the 0 2 absorption is ignored, unacceptable 
errors in the water-leaving reflectance will result (Table 5). However, given an aerosol profile — 
the base profile — through simulations we have found a simple empirical relationship (Eq. (6)) 
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that removes the effect of the absorption from the aerosol component of the reflectance (p a + p ra )* 
Combining this with a similar empirical expression (Eq. (4)) for the Rayleigh scattering component 
(p r ), allows application of the atmospheric correction algorithm developed by Gordon and Wang, 9 
that ignored the O 2 absorption. Test simulations in which the actual aerosol profile and the base 
profiles were identical above the marine boundary layer suggest that the difference in the error in 
the retrieval of the water-leaving reflectance due to the presence of the O 2 absorption is usually 
reduced to 0.001 using the empirical relations. An important feature of this method of correction 
for the effects of the 0 2 absorption is that it employs lookup tables that can be prepared in the 
absence of 0 2 absorption using more efficient radiative transfer codes. 

The dependence of the 0 2 absorption correction on an assumed profile is unsatisfying; however, 
the dependence is very weak (Figures 6, 7, and 11). One exception is the addition of volcanic 
aerosol which, if ejected into the stratosphere, will cause serious difficulty in the remo\al of the 
0 2 absorption effect (Figure 12). A second exception is the presence of thin cirrus clouds. These 
will produce effects similar to volcanic aerosols. However, in either case the atmospheric correction 
even in the absence of 0 2 absorption is seen to be degraded by the introduction of a high-altitude 
aerosol layer (e.g., compare Figures 10a, 12c, and I2d). Recalling that in the Gordon and Wang 
correction algorithm it was assumed that all of the aerosol is in the marine boundary layer, this 
suggests that the basic algorithm itself may have to be modified to incorporate the effects of such 
aerosols. A simple modification would be to use a more realistic vertical distribution of aerosol which 
would include some aerosol in the free troposphere and stratosphere. This would also improve the 
performance of the algorithm in the presence of 0 2 absorption. 

Finally, it may be possible to identify the presence of volcanic aerosol or cirrus clouds from 
SeaWiFS imagery. A method we are examining involves the SeaWiFS band at 670 nm. For C < 1 
mg/m 3 , p w at 670 nm is very small, i.e., < 0.0008. Simulations at 670 nm with volcanic aerosol 
(similar to those presented in Figure 12) suggest that the recovered p w there wuU be negative if there 
is a significant volcanic aerosol layer, e.g., a stratospheric aerosol optical thickness of 0.05 — 0.10. 
Thin cirrus clouds w T ould produce a similar effect. 
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Figure Captions 


Figure 1. A portion of the 0 2 A band at 762 run covering approximately 1.165 run. 

Figure 2. Schematic of the aerosol profile for the atmospheric model described in the text. 

I 

Figure 3. Fractional change in reflectance ( p T — p' r )/ p' T , as a function of M for a spectral band 
from 745 to 785 nm. 

Figure 5. Fractional change in reflectance (p A - p' A )/p' A , as a function of M for a spectral 
band from 745 to 785 nm. 

Figure 6. Error in the prediction of p A from p' A when the actual aerosol profile differs from 
the base profile, h is the aerosol scale height. All test profiles use M50. 

Figure 7. Error in the prediction of p A from p 1 ^ when the actual aerosol profile differs from 
the base profile, h is the aerosol scale height. All test profiles use T50. 

Figure 8. (p A - p' A )/p' A as a function of M for the selected base profile. 

Figure 9a. Error in the retrieved t(443)p lu (443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and V R = 40 km. 

Figure 9b. Error in the retrieved t(443)p UJ (443) for viewing at the center of the scan with 
aerosol profile T80-T50-B and V R = 40 km. 

Figure 9c. Error in the retrieved t(443)p lu (443) for viewing at the edge of the scan with aerosol 
profile M80-T50-B and V R = 40 km. 

Figure 9d. Error in the retrieved t(443)/>u,(443) for viewing at the edge of the scan with aerosol 
profile T80-T50-B and VR = 40 km. 

Figure 10a. Error in the retrieved t(443)pu,(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and V R — 15 km. 
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Figure 10b. Error in the retrieved £(443)/?u,(443) for viewing at the center of the scan with 
aerosol profile T80-T50-B and VR = 15 km. 

Figure 10c. Error in the retrieved £(443)p w (443) for viewing at the edge of the scan with 
aerosol profile M80-T50-B and VR = 15 km. 

Figure lOd. Error in the retrieved f(443)/> w (443) for viewing at the edge of the scan with 
aerosol profile T80-T50-B and VR = 15 km. 

Figure 11a. Error in the retrieved £(443)^.(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and VR — 40 km, and a simplified base profile. 

Figure lib. Error in the retrieved £(443)^,(443) for viewing at the center of the scan with 
aerosol profile T80-T50-B and VR = 40 km, and a simplified base profile. 

Figure 11c. Error in the retrieved £(443)^,(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and VR = 15 km, and a simplified base profile. 

Figure lid. Error in the retrieved £(443)^^(443) for viewing at the center of the scan with 
aerosol profile T80-T50-B and VR = 40 km, and a simplified base profile. 

Figure 12a. Error in the retrieved £(443)p^(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and VR = 40 km, and the standard base profile, when Volcanic aerosol 
(V) is added with optical thickness 0.05 at 550 nm. 

Figure 12b. Error in the retrieved £(443)^^(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and VR = 40 km, and the standard base profile, when Volcanic aerosol 
(V) is added with optical thickness 0.10 at 550 nm. 

Figure 12c. Error in the retrieved £(443)p^(443) for viewing at the center of the scan w r ith 
aerosol profile M80-T50-B and VR = 15 km, an d the standard base profile, when Volcanic aerosol 
(V) is added with optical thickness 0.05 at 550 nm. 
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Figure 12d. Error in the retrieved t(443)pu,(443) for viewing at the center of the scan with 
aerosol profile M80-T50-B and V R = 15 km, and the standard base profile, when Volcanic aerosol 
(V) is added with optical thickness 0.10 at 550 nm. 
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Table 1: Total aerosol optical thickness for model M80-T50-B 


VR (km) 

A (nm) 


443 

550 

765 

865 

15 

0.7029 

0.6345 

0.5509 

0.5334 

•'40 

0.3594 

0.3085 

0.2462 

0.2259 


Table 2: Boundary-layer aerosol optical thickness for model M80-T50-B 


VR (km) 

A (nm) 


443 

550 

765 

865 

15 

0.5242 

0.4974 

0.4648 

0.4538 

40 

0.1806 

0.1738 

0.1602 

0.1563 


Table 3: Total aerosol optical thickness for model TS0-T50-B 


VR (km) 

A (nm) 


443 

550 

765 

865 

15 

0.8226 

0.6345 

0.4037 

0.3289 

40 

0.4006 

0.3085 

0.1S55 

0.1589 


Table 4: Boundary-layer aerosol optical thickness for model T80-T50-B 


VR (km) 

A (nm) 


443 

550 

765 

865 

15 

0.6438 

0.4974 

0.3176 

0.2592 

40 

0.2218 

0.1738 

0.1095 

0.0893 
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Table 5: A p when existence of O 2 absorption is ignored. 


Model 

VR (km) 

Scan Center 

Scan Edge 


M80-T50-B 

40 

~0.012 

~0.014 


M80-T50-B 

15 


~0.016 


T80-T50-B 

40 

-0.028 

~0.020 


T80-T50-B 

15 
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